Atomic-scale perspective on the origin of attractive step interactions on Si(113) 



Cristian V. Ciobanu, Dhananjay T. Tambe and Vivek B. Shenoy 
Division of Engineering, Brown University, Providence, RI 02912 



Ames Laboratory 



Cai-Zhuang Wang and Kai-Ming Ho 
USDOE and Department of Physics, Iowa State University, Ames IA 50011 
(Dated: February 6, 2008) 



o 
o 

(N 

> 

o 
O 

o 

m 

(N 
> 

in 

o 
m 
o 



c 

o 
o 



Recent experiments have shown that steps on Si(113) surfaces self-organize into bunches due to 
a competition between long-range repulsive and short-range attractive interactions. Using empir- 
ical and tight-binding interatomic potentials, we investigate the physical origin of the short-range 
attraction, and report the formation and interaction energies of steps. We find that the short-range 
attraction between steps is due to the annihilation of force monopoles at their edges as they combine 
to form bunches. Our results for the strengths of the attractive interactions are consistent with the 
values determined from experimental studies on kinetics of faceting. 



Self-assembly of steps is an important way to achieve 
surface patterning at length scales where the usual litho- 
graphic techniques are not applicable. On certain vici- 
nal semiconductor surfaces, the equilibrium bunching of 
steps leads to the formation of periodic hill-and-valley 
patterns (or grooves), that can serve as templates for 
growing regular arrays of one-dimensional (ID) nanoscale 
devices like quantum wires. Experimental work in recent 
years 0, 0> El 0> Et has revealed that such groove struc- 
tures are formed on vicinal Si(113), where the flat (113) 
surface coexists with step-bunched regions over a wide 
range of temperatures and misorientation angles. 

In a series of experiments, Mochrie and coworkers 
[H Q have studied the orientational phase diagrams 
for stepped Si(113) surfaces miscut towards several crys- 
tallographic orientations. In all the cases, their stud- 
ies strongly suggest that, at low temperatures, the steps 
gather in bunches as a result of attractive step-step in- 
teractions. Independent experimental work on the fluc- 
tuations of the single and multiple-height steps and the 
kinetics of faceting transitions by Sudoh and coworkers 
0, have confirmed the presence of short-range attrac- 
tive interactions between steps. Furthermore, their anal- 
ysis of the step-zipping process has provided quantita- 
tive estimates for the magnitude of the attractive inter- 
actions [f| . The presence of such interactions has been as- 
sumed in theoretical models that show the orientational 
phase diagram Q can be explained in terms of competing 
short-range attractive and long-range repulsive step in- 
teractions 0, 0] • While much progress has been made in 
understanding the self assembly of steps on Si(113) sur- 
faces, a microscopic analysis of the step interactions has 
not been attempted and the physical origin of the short- 
range attraction has not been convincingly explained. 

In this letter, we investigate the interactions between 
steps on Si(113), with the goal of understanding the ori- 
gin of the effective short-range attraction between steps. 
Describing each step as a pair of equal and opposite force 
monopoles located at the step edges, we find that the 
short-range attraction between steps is due to the can- 



cellation of force monopoles that belong to adjacent steps 
in a step bunch. This observation is confirmed by inves- 
tigating the scaling of the step formation energies and 
interaction strengths with the step height. Noting that 
the structure of the steps is determined by the direction 
of the miscut, we consider here Si(113) surfaces miscut 
towards [332], so that our estimates for the magnitudes 
of step interactions can be directly compared with the 
recent experiments of Sudoh et al. . 

The slab geometry that shows the Si(113) terraces and 
types of steps studied in the present work is given in 
Fig. 2] Steps run along the [110] direction and consist 
of (114) nanofacets 0, with a (2 x 1) reconstruction ||. 
For the (113) surface we consider both the (3x1) and 
the (3x2) structures that are stabilized by interstitial 
atoms. We distinguish two possible structures for the 
(3x2) case, depending on the position of the interstitial 
rows on the terraces with respect to the step. When the 
interstitials on the lower terrace are away from the step 
or close to the step, we call them the (3x2)a and the 
(3x2)c structures, respectively (refer to Fig.^). 

The dimension of the simulation cell in the x-direction, 
L x , is determined by the number k of Si(113) unit cells 
and the number n of Si(114) unit cells that constitute 
the step: L x = nA + ka x , where a x — 12.73A and A = 
16.2lA are, respectively, the periodic lengths of the (113) 
and the (114) surfaces along [332] (Fig.QJ. The periodic 
length in the [110] direction L y is the smallest length 
compatible with the periodicity of both the Si(113) and 
the Si(114) structures, L y = 23.04A. The thickness of the 
slab L z was chosen to be ^150 A. To simulate the vicinal 
geo metry, we use shifted periodic boundary conditions 
[10| in which the shift along [113] is determined by the 
step height nh, where h = 1.64 A is the height of a single 
step. 

The energies of the slabs in Fig. Q are computed using 
empirical and tight-binding (TB) models for the atomic 
interactions [ill ]. Empirical calculations are carried out 
with the Tersoff potential T2 ^3] as it reproduces the 
structure, surface stress and relative energy of the Si(113) 
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FIG. 1: (a) Geometry of a typical stepped Si(113) surface 
with the simulation cell enclosed between thick dashed lines. 
Different step structures with n = k = 1 (miscut 6 = 3.2°) 
are shown on the right. Arrows mark the positions of the 
interstitial atoms that stabilize the terraces Q. For the (3x2) 
surfaces, there are two possible ways of creating a step: the 
interstitial rows on the lower terrace can be located either 
away from the step, (3x2)a, or close to the step, (3x2)c. The 
symbols di and di label the dimers on the (114) step. 

and Si(114) reconstructions reasonably well. While the 
tight-binding methods are more accurate, they are usu- 
ally impractical for ledge energy calculations since they 
require slabs with several thousands of atoms to capture 
the long-range interactions between steps. To circumvent 
this problem, we have developed a multi-scale approach 
where the tight-binding model is coupled to an empirical 
potential as described in the following paragraph. 

The key idea of the coupling scheme is to handle the 
long-range elastic fields produced by steps using an em- 
pirical potential whose elastic constants are very close 
to their TB counterparts; this would ensure that the 
elastic fields in the hybrid scheme would not be very 
different from the fields obtained in the case where all 
the atoms were treated using TB. We use a charge self- 
consistent TB method, whose superior transferability has 
been studied in detail in [l^. This method is coupled 
with the Tcrsoff (T3) potential 0| using the scheme 
shown in Fig. We divide the simulation cell into 
three regions: (1) a thin surface region where the com- 
plex bonding topologies in the vicinity of the step-edge 
are treated with the TB method, (2) a thick "bulk- like" 
region where all the atoms with tetragonal bonding en- 
vironments are treated with the T3 method and (3) a 
padding zone which is used to communicate atomic dis- 
placements between the first two regions. The atoms in 
the simulation cell are relaxed using the following pro- 
cedure (refer to Fig. 0): (1) relax the thin slab (Zone 
1) Fig- Eta)) using TB, while keeping the padding zone 
fixed; (2) update the atomic coordinates of Zone 1 in 
Fig. [2£b) with the values obtained after the TB relax- 
ation of Zone 1 in Fig.|2Ja); (3) relax the thick slab and 
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FIG. 2: Relaxation scheme employed to compute ledge ener- 
gies using a coupling of the tight-binding (TB) method [l3j | 
with the T3 potential 1141 . Zone 1 (thin) consists of a few sur- 
face layers that are relaxed using TB (a), but are kept fixed 
during the relaxation of the thick slab (b). Zone 2 (thick) 
represents the deeper layers, which interact with the surface 
layers through a padding zone P.Z. The arrows show how 
atomic coordinates are updated during stages (2) and (4) of 
the procedure described in text. 

the padding zone (Zone 2 and P.Z. in Fig. Hfb)) using 
the T3 potential while keeping the atoms in Zone 1 fixed; 
(4) update the coordinates of the P.Z. atoms in Fig.|5fa) 
with the values obtained after the T3 relaxation of P.Z. 
in Fig.|3Jb); (5) repeat steps (1) through (4) until energy 
convergence is achieved |l5|. 

Once the total energy E of the slab is calculated, step 
interactions can be obtained from the ledge energy [Io| 

\{L x ) = {E-Ne b - lll3 L x L y )/L y , (1) 

where e& is the bulk cohesion energy per atom jl(| , N is 
the number of atoms in the cell in Fig.^ an d 7113 is the 
surface energy of Si(l 13) . Using the fact that the steps 
interact through long-range dipolar fields [l^. we have 
analyzed the short-range interactions by fitting the ledge 
energies in Eq. (JIJ to the form Q 

7T 2 A 2 

\{L X = nA + ka x ) = A„ + G n — — - U n 5 k0 , (2) 

where A„ is the formation energy per unit length of a 
step of height nh, G n is the strength of the long-range 
dipolar interactions, and Sko denotes the Kronecker delta 
symbol, where k is the number of (113) unit cells on the 
terrace (Fig. QJ. The energy offset U n that appears in 
Eq. (J2J describes the effective interactions at the small- 
est possible spacing between the steps, nA (note that the 
(113) terraces are absent in this case). A positive (nega- 
tive) value of U n indicates the presence of a short-range 
attraction (repulsion) between steps. Our results for dif- 
ferent terrace and step structures are discussed below. 

We first focus attention on the vicinal surfaces made of 
single-height steps and illustrate the general features of 
the step interactions by considering the energetics of the 
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FIG. 3: Ledge energy A of the single-height (114) steps sepa- 
rated by (3x2)c terraces, plotted as a function of The 
points corresponding to separations L x > A =16.2lA are fit 
to a straight line, A = Ai +7r 2 GiA 2 /6L^. At the smallest sep- 
aration A, the ledge energy falls below the straight line, indi- 
cating a short-range attraction U\ = 7r 2 Gi/6 + Ai — A(7i 14 — 
7113) [Ts|. The inset (a) shows an array of steps (simulated 
via periodic boundary conditions) where each step is repre- 
sented schematically as a pair of force monopoles F+ and F- . 
When L x = A, there are no terraces between steps and the 
monopoles corresponding to adjacent steps cancel each other, 
as shown in inset (b). 



(3x2)c steps obtained from the T2 potential. In Fig. [3] 
we plot the ledge energy as function l/L x (L x is the step 
separation) and provide the fitting parameters A i; G\ 
and Ui defined in Eq. J5J. If the dipolar description of 
step interactions is assumed to hold over the entire range 
of step separations (as implied by Eq. (01 ), Fig. shows 
the presence of a short-range attraction (Ui), since the 
point corresponding to L x = A lies below the linear fit 
for the ledge energy. The physical reason for the origin 
of this attraction will be discussed next. 

We start by noting that the Si(l 13) and Si(114) facets 
are atomically compatible, in the sense that they can 
meet to form edges parallel to [110] without any changes 
in their individual bonding topologies. Subsequently, the 
elastic fields of a (114) step (or nanofacet) can be ex- 
pressed in terms of a dipole which consists of a pair of 
equal and opposite force monopoles (F + , F-) 19] at the 
edges of the step, as shown in inset (a) of Fig. While 
such steps are expected to show inverse-square repulsive 
interactions at large separations, we find that this behav- 
ior holds for separations as small as 29A (this point is 
marked in Fig. EJ . When the steps are closer than 29 A, 
they lose their individuality forming a smooth Si(114) 
surface: at a separation of L x = A there is a cancella- 
tion of the monopoles at the step-edges as illustrated in 
inset (b) of Fig. The absence of dipoles with mutual 
repulsive interactions is the reason for the presence of an 
effective short-range attraction between steps. Further 
evidence for this argument is obtained by computing the 
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TABLE I: Interaction strengths and step formation energies 
for single-height steps computed using the multiscale tight- 
binding/Tersoff scheme (TB-T3) and the T2 potential 
The labels "in" and "out" denote the relative tilting (in- and 
out-of-phase) of the dimer rows di and 0I2 shown in Fig. 
The last column shows the quantity e = Ai — A (7^4 — 7113), 
which allows for a direct comparison with the experimental 
results (Exp. |f|). All quantities (except the ratio U1/G1) are 
given in meV/A. 



interactions of multiple-height steps. We will focus on 
such calculations after we discuss the interaction param- 
eters for other types of steps given in Fig. 

Using both the T2 potential and the TB-T3 method, 
we find the existence of short-range attractive interac- 
tions between steps for all the reconstructed surfaces con- 
sidered here. However, the magnitudes of these interac- 
tions depend strongly on the structures of the terraces 
and steps (refer to Table [IJ. Since the TB method al- 
lows for the tilting of the dimers d\ and di on the step 
(marked in Fig. 01, we have reported the results for the 
cases where the tilting patterns are in-phase and out- 
of-phase j2jj. The magnitude of the attractive interac- 
tion Ui is « 60meV/A and w lOmeV/A for the TB-T3 
method and T2 potential, respectively. This discrepancy 
between the values of U\ computed with TB-T3 and with 
T2 is due to the differences in step formation energies ob- 
tained with the two methods, as well as to the sensitivity 
of Z7j with respect to changes in (7^4 — 7113) (refer to 
Fig. gj. For both T2 and TB-T3, the repulsion strength 
G\ is of the order of lOmeV/A, while the formation en- 
ergy Ai is of the order of lOOmeV/A (see Table HJ. In 
order to compare our results with the experimental ones 
we note that Sudoh et al. Q report the quantity 
e = Ai — A (7^ 4 — 7113), i.e. the difference between the 
step formation energy Ai and the excess energy of the 
Si(114) facet per length of the step. Table Q] shows that 
the values for e obtained with the TB-T3 method are of 
the same order of magnitude as the estimates of Sudoh et 
al. |5|. Closer agreement could be achieved if the long- 
range interactions are included in the analysis of kinetics 
of step-zipping reported in Also, future experiments 
aimed at determining the structure of the step (a or c) 
could verify the trends given in Table [I] 

We will now focus on the formation energy and inter- 
actions of steps with height nh (n-bunches) for n = 2, 3. 
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TABLE II: Step-height dependence of the formation energy 
An and interaction parameters G n , U n for the (3x2)c struc- 
ture, determined with the Tersoff potential T2 |12|. All quan- 
tities are given in meV/A. 

Because of the larger size of the simulation cells, we have 
used only the T2 potential [12J for this analysis. In the 
zero temperature limit, the formation energy of an n- 
bunch can be expressed as the sum of the formation and 
interaction energies of the steps within the bunch Q : 

n-l 

A„=nA 1 -C/ 1 (n-l) + G 1 ^^. (3) 

i=l 

Due to the presence of attractive interactions between 
single-height steps, the formation energy of an n-height 
step is expected to show sub-linear dependence of n, 
A„ < nAi |7j. This prediction is indeed borne out by 
the atomistic results given in Table ITT1 Using the param- 
eters Ai, Gi, U\ from Table [I] for the (3x2)c structure, 
we find that the bunch energies predicted by Eq. © are 
within 1-2 meV/A of the values for A n obtained in the 
present simulations. 

Turning our attention to the step-step interactions, we 
find that the repulsive interactions G n show a nearly 
quadratic dependence of the step height nh, while the 
strength of the attractive interaction U n shows a weak 
dependence [2l| on n (refer to Table In)) . The quadratic 
behavior of G n provides strong evidence to the idea re- 
garding the cancellation of the force monopoles associ- 
ated with the adjacent steps in an n-bunch. When such 
cancellation takes place, the only remaining monopoles 
are associated with the ends of the bunch. Consequently, 
each n-bunch behaves as a force-dipole with a moment 
that is proportional to nh. Identical n-bunches interact 
as dipoles [17|, with a repulsion strength determined by 
the square of the dipolar moment of one ro-bunch, which 
results in the n 2 dependence of G n . 

In conclusion, we have shown that the effective short- 
range attraction experienced by steps on Si(l 13) can be 
explained by the cancellation of force monopoles associ- 
ated with adjacent steps in a step bunch. While numer- 
ical estimates for the interaction parameters depend on 
the interatomic potential used, our physical reasoning for 
the origin of attractive interaction is robust and is con- 
sistent with experiments. A key point in understanding 
how the atomic structure can lead to step-step attrac- 
tive interactions is the atomic compatibility of the two 
stable surface orientations, Si(113) and Si(114): when 
nanofacets with these orientations intersect forming a 



step edge, no bond breaking or rebonding occurs. Re- 
cent experimental work |2^| indicates that no rebonding 
occurs at the step-edges formed between the Si(113) and 
Si(337) facets; our calculations show the presence of at- 
tractive interactions between steps in this case, further 
lending support to our physical picture. The difference 
between the type of steps considered here, and the most 
common steps on Si(001) surfaces, such as the SB and the 
DB steps |23j is that when steps on Si(113) merge they 
form a new facet (e.g. Si(114)), while such atomically 
smooth facets are not observed in the latter case due to 
significant rebonding at the step-edges |2^|. Therefore, 
a simple cancellation of force monopoles is not expected 
for the SB and the DB steps on Si(OOl), and hence they 
do not show attractive interactions [10(. Further theo- 
retical and experimental studies aimed at elucidating the 
dependence of the step-step interactions of the azimuthal 
orientation will shed more light on the self-assembly of 
steps on Si(113) miscut in an arbitrary direction. 
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